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ABSTRACT 


The object of this thesis is an examination of the validity of 
geostrophic calculations of the downstream component of the time- 
averaged Florida Current in the Straits of Florida by comparison of 
calculated and directly measured current fields. The study is 
motivated by the assumption made in modern inertial current theory 
that downstream current speed is in geostrophic balance, plus 
evidence of recent studies indicating that this region of the 
Florida Current is primarily inertial in nature, The principal 
conclusion reached is that geostrophic calculations yield a valid 
first order approximation to the observed velocity fields, indicating 
that the assumption made in inertial current theory is valid. 

In addition, it has been shown that in geostrophic calculations in 
this region, the density field may be approximated as a parabolic 
function of temperature only. A general discussion of mass field 


adjustment to downstream speed changes is offered, 
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I. INTRODUCTION 


The purpose of this thesis is the examination of the validity 
of the geostrophic approximation for the Florida Current (within the 
Straits of Florida). An experiment was conducted during 1965-1966 
having as one objective the critical investigation of such an approxi- 
mation. Earlier studies of the application of geostrophy in this area 
have been made in order to calculate the velocity field for the pur- 
pose of examining various features of the current. Since the 
development of the free instrument technique (Richardson and Schmitz, 
1965), rapid, direct measurement of the current field is possible, and 
geostrophic calculations are no longer necessary for current field 
determination. In this thesis, interest is focused on the extent to 
which the geostrophic approximation describes the downstream component 
of the flow as postulated in contemporary inertial current theory 
(Robinson, 1965). Robinson (ibid.) further points out that use of 
Space-temperature (x,y,T) coordinates simplifies the mathematics 
involved in developing models. The geostrophic equation as examined 
is unusual in its use of an equation of state with density expressed 
as a function of temperature only. 

This thesis constitutes the first quantitative validity deter- 
mination of geostrophic calculations based on simultaneous hydrographic 
data for a steady (time-averaged) current. The validity of the 


geostrophic equation will be determined by comparing calculated 





downstream velocity fields with observed velocity fields obtained 
using the free instrument technique. These observed fields will 
represent a time-averaged Florida Current at four different sections 
across the Straits of Florida (Figure 1). 

If the geostrophic approximation is valid for the Straits, it 
will give meaningful results throughout the area of investigation. 
Direct measurement of the current shows downstream acceleration 
(Clausner, 1967). For geostrophy to hold along the length of the 
current, the mass field structure must adjust as changes occur in 
the downstream velocity fields, maintaining the balance between 
pressure gradient and coriolis force. This mass field adjustment 
will be examined in conjunction with the changing velocity fields. 

Several previous investigations have invoked geostrophy in the 
study of the Florida Current (see for example: Wiist, 1924; Parr, 
1937). The extensiveness of the free instrument data available from 
this experiment (four cross-sections over a 225 Km. downstream scale) 
plus the quality of this data afford a good basis for a critique of 
earlier investigations. Such a critique will be included as a second 


objective of this thesis. 





FIGURE 1: Chart of Section Locations. 
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A. General 

Four sections were selected at locations which provide various 
combinations of downstream scale and channel geography. It was felt 
that the geostrophic approximation could be thoroughly tested when 
applied in such a variety of situations. Furthermore, the Florida 
Straits over the length where sections were sampled is a closed 
channel, except for the Santaren and Northwest Providence Channel 
openings. This configuration permits only minimal distortion of the 
Florida Current by other sources of water. 

Section I, between Marathon (Vaca Key), Florida and Cay Sal Bank, 
B.I. has a "V" shaped bottom profile reaching depths in excess of 
1000 m. Section II is located about 100 Km. downstream between Fowey 
Rocks and Cat Cay, B.I. past the point where the axis of the stream 
has turned toward the north. At this section maximum depth has reduced 
to just over 800 m. and the bottom of the trough has flattened. 
Current width at the surface has been reduced from over 100 Km. at 
Bection I to about” &5 Km. at this section. Section Ll) is 25° Km. 
further downstream and has much the same bottom profile as Section II. 
Current width at the surface, however, has decreased slightly to about 
80 Km. Section IV is another 100 Km. downstream, from Fort Pierce, 
Florida to Matanilla Shoal, northwest of the Little Bahama Bank. The 


profile here has regained its sharp "V" profile while current width 





at the surface has diverged to 85-90 Km., and maximum depth has 
decreased further to /00 m. 

Each section consists of twelve (Sections I, IV) or thirteen 
(Sections II, III) stations. Measurements were taken at these stations 
using the free instrument technique. This technique yields data in 
the form of down and cross-stream current fields and vertical tempera- 
ture profiles (see Tables I-IV). Each of the 50 stations was sampled 
over the full surface to bottom range an average of six to twelve 
times; each sampling consisted of between one and six instrument 
drops (depending on station depth), plus a surface current measurement. 

Horizontal station spacing never exceeds 10 Km., and the spacing 
is reduced near boundaries and in the cyclonic zone where there is 
intensified velocity shear. Sampling of each section was conducted 
over a period of at least three weeks (except for Section IV, where 
only seventeen days of sampling were conducted due to equipment 
difficulties). It was hoped, by judicious sampling over this period, 
to minimize tidal biasing of the data and construct a time-averaged 
representation of the mean, steady state Florida Current at the 
selected stations. 

B. Analysis 

The free instrument method yields data in a form which can be 
applied in geostrophic calculations in a space-temperature (x,y,T) 
reference frame. The x coordinate is horizontal, has its origin at 
the left shore of the section (looking downstream) and increases 
cross-stream. The y coordinate is horizontal, parallels the axis 
of the stream, and is positive in the downstream direction, with its 


zero line along the section. The T coordinate is vertical, 





positive upward. 

Of the raw data obtained from the free instrument measurements. 
the downstream component of the time-averaged velocity (Vobs) and the 
temperature profile for each station, depths of selected isotherms 
were obtained. These have been used to construct smooth isotherm 
profiles for the four sections. Perusal of the isotherm profiles 
suggested selection of the 26° isotherm as the upper bound, as it 
stays fairly shallow (never deeper than 130 m.) and rises to the 
surface only at the extreme left hand edge of the stream. Likewise, 
the 8° isotherm was chosen as the lower bound because it is present 
across all of the sections at depths sufficient to include areas of 
relatively low velocity (less than 20 cm/sec). Intermediate isotherms 
of 10, 14, 18 and 22°C were chosen for convenience. 

Profiles of Vobs have been constructed in the same general way 
as the isotherms, but from vertical velocity profiles for individual 
stations. These contours are seen as dashed lines in all figures 
showing the geostrophically calculated isotachs (Figures 2-7). 

A detailed derivation of the form of the geostrophic approximation 
used is given in Appendix A. In brief, one starts with the geostrophic 
equation for the downstream velocity component, (v): 

O>fv = dP/ dx 
where f is the coriolos parameter (2Qsingd) and P is pressure. Q is 
the angular speed of the earth, and 6 is latitude. Taking the deri- 
vative of both sides with respect to depth, and applying the 
hydrostatic approximation: 
Pofdav/dz = gop/dx. 


Transforming to (x,y,T) coordinates and transposing, 


dv/dT = g/oof (dp/9T) (dD/dx) 
where D is isotherm depth. 

Isotherm slope dD/dx is computed for each of the six chosen 
isotherms at each station. Next the thermal shear dv/0T is plotted 
along the selected isotherms, and velocity differences between 
isotherms determined planimetrically. Starting from an isotherm 
along which downstream component of velocity is known from direct 
measurements (usually the 8°C isotherm is chosen), the downstream 
component of the geostrophic velocity field is calculated by numerical 
integration of dv/dT values. 

At this point a density-temperature relationship must be 
determined. It is desired to construct a relationship of the simplest 
form capable of producing an acceptable reproduction of the existing 
velocity field when used in the geostrophic calculations. Three 
equations of state are compared in the extent to which they achieve 


this end. One equation of state used is of the form: 


Q 
I 
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where the coefficients 


Ay = Ay (x). 
This equation will be called the variable coefficient polynomial 
equation (herein abbreviated VCP). The coefficients (Ay) have been 
determined from hydrographic data by fitting the data to a least 
squares polynomial of order K, where K is the smallest integer such 
that the RMS deviation of the fit has a maximum of .106,. The 
hydrographic data used was taken by University of Miami Marine 


Laboratory personnel at 8 stations in the Fowey Rocks - Gun Cay region, 





and has been presented in the form of a table of coefficients by 
Schmitz and Richardson (1966). Interpolation has been made where 
necessary to account for differences between the locations of the 
hydrographic and the free instrument stations. 


Another p-T relationship used is a parabloic equation: 


PeeeOutro(T-T le 


where 


one ce (a@i-T )). 
The zero subscript denotes reference quantities (i.e. based on the 
8°C isotherm). The constants og and « are determined on the basis of 
mean density characteristics of Florida Current water. 

A linear equation of state was applied in the calculations 
at one section. 

Peo ape tier lark) 

The constant y was determined as were do and kK. 
C. Errors 

All error estimates discussed in this section are considered 
upper bounds. It has been shown (Clausner, op. cit.) that errors in 
observed velocity are up to 5Z, and errors in isotherm depth are 
also up to 5%. Since the error in x is of the order of tens of meters 
and isotherm slopes are calculated over a minimum distance of 5 Km., 
errors in isotherm slope dD/0x due to errors in x are negligible. 
Isotherm slope errors are, then twice the error in isotherm depth 
or 10%. 

Errors in the d0/d3T calculation depend upon the equation of 
State. For the VCP equation, a 3-5% error is introduced; the parabolic 


equation yields errors up to about 10%, and; the linear equation gives 





10 


typical errors up to about 50%. Conditions are considerably worse at 
the left edge of the stream than in the middle or right sections, but 
this area of extreme errors is disregarded in the above estimates. 

Errors in thermal shear are merely accumulations of the dp/dT 
and 3D/dx errors. These values would be 15% for the polynomial 
equations, 20% for the parabolic equation, and 60% for the linear 
equation. Final geostrophic velocity is estimated to be in error 
by about the same amount as the thermal shear, with systematic errors 
increasing away from the reference isotherm. 

The validity of the geostrophic approximation is determined by 
calculating the mean and root mean square deviations of the geostrophic 
velocity from the observed velocity at intervals along the observed 
isotachs. The calculations are of the form: 


N 


A = 1/N 2(A,) dmg = [1L/N £(42)] 


Xs 


where A; indicates the derivation at the selected points, The 
deviations will be expressed as percentages. 

In relation to any discussion of the Florida Current as a 
geostrophic current, one should be cognizant that Webster (1962) 
has shown that horizontal gradients of Reynolds stresses produce non- 
geostrophic components of velocity which average 10% of the velocity, 


and may reach a maximum of 25% of the velocity. 
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III. RESULTS 


A. Geostrophic Approximation 

Results of the geostrophic calculations are presented as isotachs 
drawn on cross-stream bathymetric profiles. Observed isotachs are 
superimposed on the profiles as dashed lines. Calculations of percent 
deviation of calculated from observed velocities disregard the area 
of large error which occurs at the far left (looking downstream) of 
the current, 

The calculated downstream component of the velocity field has 
been computed for comparison on three different bases: first, three 
different equations of state have been used; second, calculations 
have been based on three different levels (isotherms) as reference. 
and; third, calculations were made at the four cross-stream sections. 

In comparing velocity fields calculated on the basis of the 
linear, parabolic and VCP equations, the 8°C isotherm has been used 
as reference (Figures 2-3) and Section II was chosen as the comparison 
section. The linear equation yields an extremely poor approximation 
to the observed field: in fact, calculated surface velocities are 
nearly double those observed. This equation of state applied in the 
geostrophic approximation is clearly unacceptable. Isotachs con- 
structed from calculations based on the parabolic equation of state 
tend to skew to the right, producing deviations which increase as one 
approaches the left extremity of the channel. This indicates that 


the “mean'' parabolic equation is more appropriate to waters in 
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mid-channel and to the right thereof. Mean errors using the parabolic 
are 10-20%. The VCP equation of state produces geostrophically 
calculated isotachs in generally good agreement with observed 
conditions. Again, there is a tendency for the calculated isotachs to 
shift to the right though the effect is less pronounced than in the 
parabolic calculations, One cause may be the different water mass on 
the left side of the channel. Differences between VCP calculated and 
observed isotachs are quite close, quantitatively to those between 
parabolic calculated and observed isotachs, except, as pointed out, on 
the far left of the channel, 

Again with Section II as the comparison section, isotachs were 
constructed based on the VCP equation of state using the 8°C, 18°C, 
and 26°C isotherms as reference (Figures 3-4). As expected the error 
increases as one moves away from the reference level. Largest 
percentage errors occur around the 20 cm/sec isotach when the 26°C 
isotherm is the reference. Minimum overall discrepancies exist when 
the 18°C isotherm is reference. 

Examining the validity of the geostrophically constructed v 
fields at the different sections (Figures 2, 3, 5-7), the parabolic 
and VCP equations of state were used with the 8°C isotherm as reference. 
Best agreement is found at Sections I, II and III. On the average, 
isotachs constructed on the basis of both the VCP and parabolic 
equations agree to about 10-20% with observed values. 

At Section IV, the largest discrepancy occurs, about 15-20%. 

This could be due to the influx of water of different characteristics 
through the Northwest Providence Channel (Finlen, 1966). In Sections 


II and IV the parabolic equation has given better results than the 
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FIGURE 2. a. GEOSTROPHIC ISOTACHS - Section II 
8°C reference isotherm 


Linear equation of state 


b. GEOSTROPHIC ISOTACHS -— Section II 
8°C reference isotherm 


Parabolic equation of state 
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FIGURE 3. a. GEOSTROPHIC ILSOTACHS - Section II 
8°C reference isotherm 


VCP equation of state 


b. GEOSTROPHIC ISOTACHS - Section II 
18°C reference isotherm 


VCP equation of state 
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FIGURE 4, 


GEOSTROPHIC ISOTACHS - Section IIL 
26°C reference isotherm 


VCP equation of state 
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FEGURE 5. 


a. GEOSTROPHIC ISOTACHS - Section I 
8°C reference isotherm 


Parabolic equation of state 


b. GEOSTROPHIC ISOTACHS = Section I 
8°C reference isotherm 


VCP equation of state 
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FIGURE 6, a. GEOSTROPHIC ISOTACHS —- Section III 
8°C reference isotherm 


Parabolic equation of state 


b. GEOSTROPHIG USO TACHS —sSectton TEL 
8°C reference isotherm 


VCP equation of state 
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FIGURE 7, a. GEOSTROPHIC ISOTACHS - Section IV 
8°C reference isotherm 


Parabolic equation of state 


b. GEOSTROPHIC ISOTACHS - Section IV 
8°C reference isotherm 


VCP equation of state 
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VCP equation for the bulk of the current. 
B. Mass Field Adjustments 

Mass field adjustment in response to changing downstream velocity 
components is displayed as a plot of average isotherm slope and 
average downstream speed along the isotherms at each of the sections 
(Figure 8). Average slope of each of the six selected isotherms is 
plotted, as is the average velocity along each isotherm. The average 
isotherm slope was determined graphically from isotherm profiles for 
the sections. Average velocity was computed as a number average of 
the observed downstream speeds along each isotherm. 

From Section I to Section II speeds increase from the upper 
waters (26°C region) down through the 10°C region. Isotherm slope, 
though less at the 26°C and 22°C lines, increases markedly below the 
22°C level indicating, for reasons of continuity, deepening of the 
isotachs as the isotherms rise. From Sections [II to III little 
difference appears in the plots, except for the increased slope of the 
10°C isotherm at Section III, which reflects moderate speed increase 
in the mid layers (22°C - 10°C). Between Sections III and IV the 
velocity profile changes quite differently. Speeds are lower at 
Section IV down to the vicinity of the 14°C isotherm, and below this, 
Section IV speeds are higher than those of Section III. Down to just 
above the 14°C isotherm, Section IV isotherm slopes are greater. 
Below this point, the Section III isotherm slopes are greater, 
especially the 10°C isotherm. This reflects the first marked 


increase in speeds along the lower isotherms. 





FIGURE 8. 


Mass field -— velocity adjustment 


Sections I-IV 
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C. Earlier Work 

The geostrophic calculations of Wltst (op. cit.) and the agreement 
of his results with the direct measurements of Pillsbury (1890) have 
done much to convince oceanographers of the value of the geostrophic 
approximation, It is remarkable that, using temperature measurements 
taken in 1878-81, some salinity measurements taken in 1914 and a T-S 
correlation based on North Atlantic waters (to cover zones where no 
salinity values were available) WUst was able to construct a density 
field which yielded such good results when compared with the Pillsbury 
measurements of 1885-1886. His results agree (in the vicinity of our 
Section II) with those obtained using the free instrument method to 
roughly 20%. It is interesting that the isotachs resulting from his 
calculations are skewed to the right (as compared with the observed 
isotachs) much as are those calculated in this thesis. 

In his analysis of five hydrographic stations between Miami and 
Bimini, B.I., Parr (op. cit.) showed the variations in the density 
field with time. In particular he described the cross-stream 
oscillations of a high salinity core. Broida (1966) has shown how 
this variation in the mass field can result in a distorted picture 
of the current field as calculated geostrophically. One example of 
such a distortion is the appearance of a bi- or multi-axial surface 
current profile cross-stream, a phenomenon which is not evident using 


the more nearly synoptic free instrument technique. 





41 


IV. SUMMARY AND CONCLUSIONS 


During the summer months of 1965 and 1966, four sections across 
the Florida Straits:‘(over a 225 Km. downstream distance) were sampled 
using the free instrument technique, over time and space scales 
designed to yield a picture of a time-averaged, steady state Florida 
Current. From this series of measurements profiles were constructed 
for each section of the downstream component of velocity and of depth 
of selected isotherms to provide a basis for evaluation of geostrophic 
velocity fields. 

Based on hydrographic data obtained in 1962-63 in the Fowey Rocks - 
Gun Cay area, three equations of state were developed in which density 
was expressed as a function of temperature only. A fifth order 
polynomial equation was used with coefficients a function of cross- 
stream distance; a simple parabolic equation was applied, and; a linear 
equation was employed at one section. Each of these equations was 
used in the geostrophic equation for downstream velocity in a space- 
temperature coordinate system, to determine the validity of this form 
of the geostrophic equation. Three different levels were tested as 
reference for the calculations. These were the 8°C, 18°C and 26°C 
isotherms, along which velocity was accurately known from the direct 
measurements. Velocity fields from the geostrophic calculations 
have been compared with those obtained by direct measurement. 


Lastly, an attempt has been made to show the adjustment of the 
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mass field (specifically the temperature field) to changes in the 
downstream component of the velocity field as the current is accelerated 
downstream, 

Results of the experiment discussed in this thesis point to the 
conclusion that, under certain conditions, the geostrophic approxi- 
mation yields a model of the downstream component of the velocity 
field which is a first order approximation to the directly observed 
field, For most of the current, geostrophically calculated velocities 
agree within experimental error with directly measured velocities, 
Necessary conditions for such a conclusion include the use of an 
accurately known reference level and use of an equation of state in 
which density is expressed in terms of temperature by a second (or 
higher) order equation of state. Constants for the equation must be 
determined on the basis of local hydrographic conditions. These 
results confirm the validity of using the geostrophic approximation 
in theoretical models of inertial flow. 

Furthermore, it has been shown that there is a direct relation- 
ship between the mass distribution (specifically the temperature) 


and the downstream component of the velocity field, 
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APPENDIX A: 


DERIVING THE PARTICULAR FORM OF 


THE GEOSTROPHIC EQUATION 


45 





46 


Robinson (ibid.) has stated that the downstream component of an 

inertial flow is in geostrophic balance. In (x,y,z) coordinates’ 
Pofv = dP/ dx 
where the use of a constant 9, is an approximation resulting in less 
than 1% error, Taking the derivative with respect to depth of the 
above equation and applying the hydrostatic approximation: 
O>of(dv/oz) = -g(dp/dx). 
One may transform to (x,y,T) coordinates by: 
(dv/ 02) x y = (dv/9T) ~ y(0T/02z)y y 
and; 
Cop /ax)y , x ~ (8p /8T) (8D/dx), pCdT/ dz), 
which yields: 
dv/dT = g/fp,(dp/9T) (oD/ dx) 

where D is the depth of an isotherm T. 

It now remains to take the temperature derivative of the 
three equations of state, which is straightforward and gives the 
results shown below. 

Linear equation; 

dp /98T = Poy 
Parabolic equation; 
op/3T = P44, [1+2« (T-To) ] 
VCP equation; 


go /aT = © [KAT OK-D 
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